Tunable front interaction and localization of periodically forced 

waves 



Catherine Crawford* and Hermann Riecke 

Department of Engineering Sciences and Applied Mathematics, Northwestern University, 

Evanston, IL 60208, USA 
(February 8, 2008) 

Abstract 



In systems that exhibit a bistabihty between nonhnear travehng waves 
and the basic state, pairs of fronts connecting these two states can form lo- 
cahzed wave pulses whose stability depends on the interaction between the 
fronts. We investigate wave pulses within the framework of coupled Ginzburg- 
Landau equations describing the traveling- wave amplitudes. We find that the 
introduction of resonant temporal forcing results in a new, tunable mechanism 
for stabilizing such wave pulses. In contrast to other localization mechanisms 
the temporal forcing can achieve localization by a repulsive as well as by an 
attractive interaction between the fronts. Systems for which the results are 
expected to be relevant include binary-mixture convection and electroconvec- 
tion in nematic liquid crystals. 
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I. INTRODUCTION 



Localized structures have been observed in a range of pattern-forming nonequilibrium 
systems. One type of localized structure occurs when one pattern is embedded within 
another pattern. Examples include the coexisting stationary domains of long and short 
wavelengths observed in Taylor- Couette flow between co-rotating cylinders in Rayleigh- 
Benard convection in narrow slots 0, and in parametrically excited waves in ferrofiuids 
0]. Solitary waves drifting through a stationary pattern are found associated with a parity- 
breaking bifurcation in directional solidification the printer instability viscous fin- 
gering 0, cellular flames [0, and Taylor- vortex flow |^. 

In this paper, we investigate a class of localized states in which the pattern is confined to a 
small region that is surrounded by the unpatterned state, or vice versa. For example, solitary 
standing waves ('oscillons') have been observed in vertically vibrated granular layers and 
colloidal suspensions ||10i . Localized traveling waves have been observed as one-dimensional 



pulses in binary-fluid mixtures [p!T|-[T5| and as two-dimensional 'worms' in electroconvection 
of nematic liquid crystals |jl6[ . 



For a general understanding of such structures the mechanisms that are responsible for 
their localization are of particular interest. A number of different types of localization 
mechanisms have been identified (e.g. [0]). To provide a context for our results we briefly 
review the main mechanisms. 

The stable coexistence of domains of long and short wavelengths can be understood to 
be due to the instability of the constant wavenumber state combined with the conservation 
of the phase [|l^-|22|. Localized patterns can also be stabilized by a non-adiabatic pinning 
of the large-scale envelope to the underlying small-scale pattern p3| , |24| . This pinning has 



been suggested as a possible localization mechanism for oscillons ||25|| . 

To understand the traveling-wave pulses in binary mixture convection, two mechanisms 
have been put forward, dispersion p6|-p8| and the advection of a slowly decaying concentra- 
tion mode [p9|-|32| . Within the framework of the complex Ginzburg-Landau equation with 
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strong dispersion, pulses and holes can be viewed as perturbed bright and dark solitons of 



the nonlinear Schrodinger equation [p3|-p5|] . For weak dispersion, pulses have been described 
as a pair of bound fronts. In the absence of dispersion they are unstable. However, disper- 
sion may result in a repulsive interaction between the two fronts and a stable pulse can arise 
P7| , |25[] . Similarly, the advected mode modifies the interaction between fronts and can also 
provide a stabilizing repulsive interaction. A similar advective mechanism has been invoked 
|36| to explain the two-dimensional localized waves ('worms') that have been observed in 
electroconvection in nematic liquid crystals . 



More generally, the coupling of a pattern to an additional undamped (or weakly damped) 



mode can lead to its localization [^]. In the drift waves arising from a parity-breaking 
bifurcation, the local wavenumber of the underlying pattern plays the role of the additional 
mode p8| , p9| . For the oscillons in vibrated granular media it has been suggested that a 



coupling of the surface wave to a mode representing the local height of the granular layer is 



important [Q. 

For traveling waves it is well known that the external application of a resonant tem- 
poral forcing excites the counter propagating wave pT|-^]. A natural question is therefore 



whether the counterpropagating wave can play a role similar to the various additional modes 
mentioned above and can thus lead to the localization of the traveling wave into a pulse. 
Since the temporal forcing is easily controlled externally this localization mechanism would 
be tunable. 

In this paper we investigate the effect of time-periodic forcing on spatially localized 
waves that arise in systems exhibiting a subcritical bifurcation to traveling waves as is, for 
instance, the case in binary-mixture convection. We expect the results also to be relevant 
for the worms observed in electroconvection in nematic liquid crystals. 

We first consider the effect of forcing on the interaction of fronts in the absence of other 
localization mechanisms and show that forcing alone can lead to localized structures. This 
new localization mechanism can stabilize pulses with either a repulsive or an attractive 
interaction. While the interaction strength between fronts is usually determined by the 



system parameters, forcing allows the strength to be controlled externally. To focus on the 
interaction of the fronts due to the temporal forcing, we start in Section |I| with two coupled 
dispersionless Ginzburg-Landau equations and derive evolution equations for the fronts. 
In Section |IIT| we discuss the resulting front equations and compare them with numerical 
calculations. Section |IV| extends the analysis and discussion to include holes and multiple 
pulses. The combined effect of temporal forcing and dispersion is investigated in Section 
and conclusions are presented in Section ^ 



II. DERIVATION OF THE FRONT EQUATIONS 



Motivated by the pulses observed in binary mixture convection |Tl|-|T5| and the worms 
in electroconvection |]TB[ we consider a subcritical bifurcation to traveling waves in a one- 
dimensional system that is parametrically forced. To obtain a weakly nonlinear description, 
physical quantities like the temperature T of the fluid in the midplane in convection, say, 
are expanded in terms of the amplitudes A and B, of the left and right traveling waves, 

r = ^^A{x,t)e'^''^^-^'^ + B{x,t)e-'^^^^+'^'^} + c.c. + h.o.t, e < 1, (1) 

where t and x are the fast time and space coordinates. The amplitudes A and B are allowed 
to vary on a slow time scale t and a slow spatial scale x. Due to the external periodic forcing 
g2 ^g«t<Jet close to twice the Hopf frequency, Ue = 2{uh + e^fi), the expansion (|l]) is chosen in 
terms of the forcing frequency. The forcing excites the oppositely traveling waves and breaks 
the time translation symmetry, resulting to lowest order in a linear coupling between the 
two wave amplitudes. Using the remaining spatial translation and reflection symmetries of 



the system, the form of the amplitude equations for A and B can be derived ||41|-^. Hence 
we study the following set of coupled Ginzburg-Landau equations as a model describing the 

system, 

dtA = -sd^A + e^d2d^^A + nA + c\A\'^A - plAl'^A 

-g\B\'^A - r\A\'^\B\'^A - u\B\^A + uB* (2) 



dtB = +sd^B + e^d2d^^B + fiB + c\B\^B - p\B\^B 

-g\A\^B - r\A\^\B\^B - u\A\^B + vA\ (3) 

where the forcing coefficient v and the group velocity s are reaL All other coefficients may 
be complex. However, to focus on the interaction of the fronts due to temporal forcing, the 
front equations are derived for the case in which all of the coefficients are real, i.e. neglecting 
dispersion and detuning. Nonlinear gradient terms which would also appear in equations 
(§) and (D) to the order considered have also been neglected. In order for the bifurcation to 
be weakly subcritical, the cubic coefficients must be small enough to allow a balance with 
the quintic terms. 

The bifurcation to spatially extended traveling waves in electroconvection in nematic 



liquid crystals appears to be supercritical |44|. To explain the observation of worms already 



below the onset of spatially extended waves, it has been argued that a weakly damped mode 
is relevant that is advected by the wave |^| . Already below threshold such a mode can lead 



to the existence of infinitely long worms, i.e. to convection structures that are narrow in 
the y-direction, say, and spatially periodic in the x-direction [p6| , p!7| . They are expected to 



arise in a secondary bifurcation off the spatially extended traveling waves |^ . Focusing on 
the dynamics of the worms in the x-direction, the head and the tail of the worms can be 
considered as leading and trailing fronts that connect the nonlinear state, which is strongly 
localized in the ?/-direction, with the basic non-convective state. It is reasonable to expect 
that Ginzburg-Landau equations for a subcritical bifurcation will capture the qualitative 
aspects of these structures. 

With the usual scaling x = ex complex Ginzburg-Landau equations are obtained in which 
the growth term ^A is balanced by the diffusive term ddxxA. For group velocities of order 
1 this implies that the term sd^A is inconsistent with the rest of the equation, i.e. it is of 
lower order. However, by considering slower spatial scales x = e^x, the advective term sdxA 
is of the same order as fiA and the diffusive term then appears in the re-scaled equations 
(Q) and (0) as a higher-order correction of C(e^) as indicated in (Q) and (H 
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We are interested in localized solutions made up of two bound fronts connecting the basic 
state with the nonlinear state as sketched in Figure ^ One contribution to the interaction 
between the fronts arises from their overlap. For large distances it is small, inducing a weak 
attractive interaction, which destabilizes pulses. However, since the diffusion is weak ((9(e^)), 
resulting in steep fronts in A, the overlap between the fronts is of higher order in e and the 
associated interaction can be ignored. Hence the interaction between the fronts is dominated 
by the presence of the oppositely traveling wave excited by the periodic forcing. The small 
diffusion coefficient causes internal layers to arise at the positions of the fronts. The two 
bound fronts are then divided into the five regions sketched in Figure |I|, where xl and xr give 
the positions of the left and right fronts. Within regions I, III, and V the amplitude of A is 
constant. Regions II and IV are regions of rapid change where the dynamics of the fronts are 
determined. The internal layers have a width Ax = 0{e), implying Ax = e~^Ax = 0{e^^), 
which is still large relative to the critical wavelength of the traveling waves. 

We consider the forcing to be small, u = 0{e), so that the corresponding amplitude 
B of the left-traveling wave that is excited by the forcing is of the same order in e. The 
amplitudes and the parameters /i and z/ are expanded as 

A = Ao + eAi + ---, B = eBi + ---, (4) 
= /ic + e^/i2 H , u = eui-l . (5) 

where Hc = —Sc^/lGp is the value of the control parameter at which a single, non-interacting 
front is stationary. For fi > fic the nonlinear convective state invades the basic state. In 
the following derivation we go into a reference frame moving with the group velocity s. The 
positions xr and xl of the right and left fronts evolve then on a slow time scale, T = e^t. 

Inserting the expansions and into equations (0) and @ we obtain at leading order 
equations for A and B in the outer regions I, HI, and V, 

= ^icAo + cAl - pAl (6) 
= 2sd^,B, + fi,B, + u,Ao - gAlB, - uA^B^. (7) 
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Solving equation (|]) results in = in regions I and V and = = 3c/ Ap in region 
III. The corresponding solution to (0) is 



Bi{x) = -Ai + K'e^- 
a 



(8) 



where j corresponds to the regions I, III, and V and = —jic + Q^AqY + u{AqY. 

In the inner regions II and IV, the solutions vary on a fast space scale x/t which is 
captured by introducing the inner coordinates rj = (x—xl) / e and r] = (xr—x) / e, respectively. 
The spatial derivative then transforms as —>■ ±"9,,/ e. The resulting leading-order equations 
for A and B are 



= d2dr,rjAo + ficAo + cA^ - pAl, 
= ±2sd^rjBi. 



(9) 
(10) 



From d^) one obtains the left and right front solution 



1 



1 + tanh f I 



(11) 



where A^ and rj are defined above. Equation (^) implies that Bi does not depend on the fast 
space variable r], i.e. Bi = Bl'^^ is constant in regions II and IV. At 0(e), no inhomogeneity 
arises in the equation for Ai, which is therefore taken to be identically 0. But at 0{t^) the 
following equation for A2 is obtained: 

^^L,Rr~, ^ , „ d2 /I , „ /l3 d2 



CA2 = T^FdvAo - /i2^o + gBfAo + rA'.Bf - i^,B, 



(12) 



where £ = (i2C?»?r) + /^c + 3cAq — bpA^ is the linearized operator. It is singular and has the 
zero-eigenmode S^Aq, which leads to a solvability condition for equation (p!2D. The result 



of projecting equation (0) onto the zero-eigenmode determines the velocity of the fronts in 
regions II and IV, 



± 



dx 



L,R 



dT 



-2/i, + {2g + Alr)B^ 



II.IV' 



2 



II JV 



Ar 



(13) 



with B^^'^^ yet undetermined. 



Since B is generated by A, we consider the case where B = ahead of the pulse (i.e. 
for a right-travehng pulse region V where A = 0). This implies that = from equation 
Now matching the inner and outer solutions for B at the left and right positions {xl 
and xr) of the fronts, one obtains the constant values B^^ = {uiAc/a) [1 — exp {—aL/2s)] 
and B^^ = 0. The constants and K^^^ are nonzero and we see from equation that 
B grows spatially to the left in region III and approaches the limiting value uiAc/a. Then, 
in region I where the pulse has passed and A again equals 0, B decays exponentially to 
zero. The individual front velocities are given by substituting the values of B^^'^^ into the 
expressions given in (|13D 



dXL 

dxR 
dT 



-2(/i - /ij + 



u^{2g + Alr)Al _ Auf 



(14) 
(15) 



Combining these results yields the following equations describing the evolution of the pulse 
length L = xr — Xl and the velocity of the pulse relative to a frame moving with the group 
velocity s in terms of a "center of mass" coordinate M = {xr + xl)/2, 



dL 

dT 
dM 



fci(/i - /^c) + (1 - e 
_ 1 
dT~ 2 



-A;2Z/^(1 - e"^^) + A;3Z/^(1 - e 2 



(16) 
(17) 



where ki = 4^, k2 = 4^/q;, and = [{2g + Alr)AlC,]/a'^. While both ki and ^2 are always 
positive, ^3 may be either positive or negative. As will be seen in the next section, the sign 
of fcs determines whether stable pulse solutions may exist. The interaction length is given 
by 2s/a. 



III. DISCUSSION OF FRONT EQUATIONS 

The possible pulse solutions of (|T6|) and (|T^ are more easily discussed in terms of the 
quantity A = [1 — exp (— aL/2s)] which monotonically increases from zero to one as L 
increases from zero to infinity. Equation (^) is then given by 



^ ■ ^ = /(A) = - /X,) + k,v'K - (18) 

Figure || depicts /(A) for the two cases jj, > Hc and /i < Hc- For /^a > the parabola opens 
downward and the maximum occurs at A = k2/2k^, independent of forcing strength v. Fixed 
points Ao are indicated where /(A) = 0: 



k2± \ kl + *^n)k.k, 
Ao = (1 - e'^^°) = ^ (19) 

Pulses with A = Aq are linearly stable if /'(Ao) < 0. If k^ < 0, the parabola is opening 
upward with the minimum at A = k2/2k^ < so that /'(A) > for all values of A > 0. 
Hence for a stable pulse to exist it is necessary that > and in this case corresponds to 
the upper branch of solutions (|T9l). More precisely, stable pulses exist as long as k^, > k2/2 



with their length diverging for k^ k2/2. Unless stated otherwise the discussion of the 
front interaction will focus on the regime < k2/2k^ < 1 where a stable branch exists. 

From the expression (^) and the parabola in Figure 0, one can easily see the effect of 
forcing on the stable pulse length. As the forcing u increases, /(A = 0) = /(A = k2/k^) = 
ki{fi — fic) remain fixed while the parabola steepens (see Figure |^). Thus for fi > fic the 
stable pulse length is larger than Lc corresponding to Ac = ^2/^3, but approaches Lc as 
forcing increases (Figure ^). Conversely, for fi < fic the pulse length increases to Lc with 
increasing z/ (Figure Qb). As forcing u decreases for fi > /ic, the pulse length increases and 
eventually diverges to infinity. In contrast, for n < Hc decreasing forcing results in the pulse 
length decreasing until the pulse is destroyed in a saddle-node bifurcation. 

Figure ^ shows both the analytical results for the steady state solutions of equation (p^Bf) 
and numerical results obtained by integrating the full amplitude equations (H) and (|^). A 
linearized Crank-Nicholson scheme was used to solve the coupled equations and the pulse 
lengths were measured at half-amplitude. The pulse length is plotted as a function of forcing 
strength for several values of the control parameter /i confirming the expected behavior. The 
top figure is for group velocity s = 20.0 and the bottom figure for s = 1.7. Here we see that 
although the analytical result no longer agrees quantitatively for smaller s, it still describes 
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the qualitative behavior of the pulse length as the control parameters are varied. This 
dependence on the group velocity will be discussed at the end of this section. 

A numerical control technique has been employed to obtain the unstable pulse solutions 
indicated by the long dashed curves in Figure ^ Since the analysis leads to a single ordinary 
differential equation (|T^) describing the evolution of a pulse, the dynamics of the pulse length 
are essentially restricted to a one- dimensional manifold. The control technique is therefore 
relatively straight-forward. Pulse lengths are measured after evolving equations (0) and 
over a short time interval and compared with the measurements taken at the previous time 
to obtain the direction and rate of growth. This information is then used along with the 
desired pulse length to adjust the control parameter accordingly. This process is repeated 
until a steady-state pulse with the specified length is obtained. In Figure ^ the parameter 
adjusted to control the length is the forcing strength v while all other parameters remain 
fixed. The same technique is also used for varying the control parameter /i (cf. Figures §b, 



I and |I2]) 



The different regimes can be understood by looking at the individual interaction terms 
in equation (|T^) and their origins from equation (|^). The first term in ( [T^ ) is a measure 
of how far the control parameter /i is from the critical value [i^ where, in the absence of 
forcing, an isolated front is stationary. It provides a "pressure" which is directed outward 
for fi > fic and which has to be balanced by the interaction terms. The second and third 
terms describe the interaction between fronts due to forcing. The second term arises from 
the linear coupling between A and B introduced by the forcing in equation (0) through 
which B excites A. It therefore enhances the invasion of the nonlinear state into the linear 
state, which corresponds to a repulsive interaction between the leading and the trailing front. 
The third term stems from the nonlinear coupling which, when > 0, suppresses A and 
therefore weakens the invasion, implying attraction between the fronts. 

The effect of forcing on a pulse depends on the distance between the fronts. Since the 
pulse is traveling to the right, the amplitude B is growing spatially to the left. For short 
pulses, therefore, B remains small and the linear coupling term dominates the interaction 
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implying a repulsive interaction. Thus, for fi < fic the inward "pressure" can be balanced by 
a repulsive interaction if the pulses are sufficiently short. With increasing pulse lengths, B 
reaches larger values at the trailing front and the nonlinear coupling term gains importance. 
Thus, in contrast to many other localization mechanisms, the forcing can induce an attractive 
interaction that grows with distance. It is able to balance the outward "pressure" for /i > /i^. 
In this regime the pulses become shorter with increased forcing. Figure § shows two stable 
pulse solutions obtained by numerically integrating equations (0) and (H). The control 
parameter /x = —1.240 for the longer pulse and /i = —1.250 for the shorter pulse. All other 
parameters are the same for both pulses. Note that the forcing u is the same for both pulses, 
yet the amplitude B has not saturated for the shorter pulse. 

From equations (|^) and (|1^) we see that when dL/dT = for a steady pulse Lq, then 
dM/dT = ki{fi — fic)/2. Thus the pulse velocity in the moving frame is given by the invasion 
speed and depends on fx. By contrast, in dispersively stable pulses the velocity depends on 
the control parameter only through the nonlinear gradient terms [^,^. Strikingly, in the 
asymptotic calculation leading to ( PTD the pulse velocity does not depend on the forcing u. 
Figure ^ shows that the linear dependence on /i is in good agreement with this analytical 
result, but that the velocity does also depend on forcing z/. We suggest that the discrepancy 
may be explained by the effect of the excited counter-propagating wave on the leading front 
at xr. Since in the weakly nonlinear regime diffusion is much weaker than advection (for 
group velocities of 0{1)), B^ = in region IV, which implies that forcing has no affect 



on the velocity of the leading front (equation |T^). However, non-zero B within this inner 
region will slightly alter the velocity of the leading front and consequently that of the pulse 
as a function of forcing. This will also impact the length of the pulse. For the parameters 
used in Figure ^, note that B varies slowly in region II whereas it varies almost as fast as 
A in region IV suggesting that the missing contribution is more prevalent at the leading 
front than at the trailing one. Figure shows single front velocities for both the leading 
and trailing front. The curves are the analytical results obtained from (|r^) for single fronts 
{L — s> oo) and the symbols indicate numerical results. As expected, the velocity of the 
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trailing front is well described by our analysis, but the velocity of the leading front does 
depend slightly on the forcing. The leading-front velocity then dictates the velocity of the 
pulse. For larger group velocity s, the interaction length 2s/ a increases and the separation 
of the fast and slow spatial scales at the positions of the fronts becomes more distinct, as 
shown in Figure ^. Thus we expect and the numerical results confirm (Figure |^) that the 
agreement with analytical results (|16|) and (|17D for pulse length and velocity improves with 
increasing group velocity s. For small forcing u the relative error of the velocity is on the 
order of the numerical accuracy. We note that the qualitative behavior of the pulse length 
is consistent even for smaller group velocity. 



IV. HOLES AND MULTIPLE PULSES 

The analysis from section |1| can be applied to other front configurations such as hole 
states and multiple pulses. A hole consists of a localized region in which the amplitude of 
the pattern is very small and which is surrounded by the nonzero traveling-wave amplitude. 
In contrast to the hole- type solutions found in the Ginzburg-Landau equation p5| , ^ , which 
are qualitatively different objects than pulses, the holes to be discussed here are similar to 
pulses in that they are also compound objects made of two fronts connecting the basic state 
with the nonlinear state. As with the pulses, the central equation for their description is an 
evolution equation for the length Lh of the hole. It is given by 

^ = -hif, - /ij - k^v\l + e-^^'O + k,v\\ + e-'-^^% (20) 

where ki,k2, and k^ are defined as before and a stable hole may be possible when k^ > 0. 
Note that the interaction length for holes is now given by 2s/\^c\ which is longer than the 
interaction length 2s/ a for pulses. Solving dLh/dT = again yields two possible solutions 
L^, with the longer hole unstable and the shorter hole stable. In contrast with the pulses, 
stable holes thus exist only over a finite range of shorter lengths. The RHS is quadratic 
in Afi = exp(— |/ic|-^^/i/2s), describing a parabola opening upward. Although varying forcing 
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now results not only in stretching the parabola, but also in shifting it vertically, the variations 
in hole length again depends on the sign of the first term in (pOl). When fi < fi^, increasing 
forcing causes stable holes to get shorter, while for fi > Hc the stable hole length will increase. 
Although the tendency of a hole to grow or shrink with increased forcing depends only on 
the sign of /i — fic, the limiting behavior is different depending on the relative sizes of k2 and 
^3. If ks < k2 then holes exist only for ^ < and as u increases, the hole length goes to 
zero. If /cs > k2/2{\^ — 1) ^ 1.2A;2, holes exist only for /i > /ic and as u increases the stable 
hole length grows and eventually disappears in a saddle-node bifurcation with the unstable 
hole. Finally, for the intermediate values of fcs, holes exist for values of fi both above and 
below /ic. In this case, as u ^ oo, the hole length approaches the limiting length determined 
when fi = fie- 

Figure ^,b shows a stable numerical hole solution and the control parameter as a 
function of the hole length for different values of the forcing. Stable and unstable solutions 
are indicated by solid and open symbols, respectively. Again, the unstable holes are ob- 
tained by means of a numerical control technique. First we note that according to (^), the 
minimum of these curves should all be at the same value of L^, but our numerical results 
show the minimum shifted to the right as u increases. The analysis requires that A vanish 
in the hole region, but the presence of B in this region actually generates small nonzero 
A, which in turn generates B. Hence, the actual value of B at the trailing front is greater 
than predicted. This suppresses A there and the trailing front slows down, leading to longer 
pulses and a shift of the minimum of the curves to the right as forcing u increases. If the 
parameters are such that the control parameter /i can be taken smaller, the basic state is 
more strongly damped. Hence, within the hole region A is smaller and the shift of the curves 
is less pronounced as seen in Figure ^ Here c = 2.58 so that fi^ ~ —1.248 as compared to 
c = 1.8 in Figure || where yUc ~ —0.608. 

Arrays of multiple fronts can be combined to form multiple pulses. For a two-pulse 
configuration, there are four inner front regions which match to five outer regions resulting 
in evolution equations for the distance L2 between the pulses as well as the widths Li and 
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Ls of the leading and the traihng pulse, respectively. This is to be contrasted with the 
description of multi-pulse solutions in the strongly dispersive case (without forcing), where 
the pulse widths can be adiabiatically eliminated in favor of the the distance and the phase 



difference between the two individual pulses [^]. The three evolution equations for -^12,3 
are given by 

^ = /ci(/i-/i,) + fc2^^'Ai -fcsz/'Ai, (21) 
^ = -fci(/i - /.,) - fc2^%(l + A2) + hu'A\{l + A^), (22) 
^ = fci(/x - /i,) + W {A1A2 + [1 + A3(AiA2 - 1)]} 

-hu'' [k\KI + [1 + A3(AiA2 - 1)]'} , (23) 

where the Aj are defined in terms of Lj as follows: Ai = 1 — exp(— aLi/2s), A2 = 
exp(/icL2/2s), and A3 = exp(— 0^3/25). Since equation ( pT]) for Li is the same as equation 
(p!6|) which describes a single pulse, the leading pulse length is unaffected by the trailing 
pulse. But the trailing pulse length depends on both the length of the leading pulse and 
the distance between the pulses and is typically shorter than the leading pulse. Solving 
equations (pT]) and ( pB]) for the fixed points Aio, A20, A30, it is found that Aio ■ A20 = ^2/^3 
and A30 = A;3 ■ (1 — h.iQ/{k^ — k2). Since Li increases with increasing Ai, and L2 and L3 
increase with decreasing A2 and A3, this suggests that all three lengths will either increase 
or decrease as the control parameters and v are varied. The eigenvalues obtained from 



linearizing equations (|2T|) and (^) indicate that ai = a^, and a two-pulse solution may be 
stable when ^3 is positive. A stable two-pulse solution is therefore expected to exist and 
be stable whenever a single pulse exists. However, since the trailing pulse is narrower than 
the leading pulse, it may become too short and collapse - as it merges in a saddle-node 
bifurcation with the shorter unstable pulse - for parameter values where the single pulse is 
still stable. Numerically, we observe that as the control parameter /i is decreased all three 



lengths get shorter and the trailing pulse eventually collapses to zero. Figure [1^ shows a sta- 
ble two-pulse solution superimposed over a single pulse solution confirming that the leading 
pulse length is unaffected by the trailing pulse. 
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V. DISPERSION EFFECTS 



Waves generally have both linear and nonlinear dispersion. Hence, the coefficients and 
the amplitudes in equations (0) and (|^) are in general complex. In [^,^, it is shown that 
for weak dispersion the interaction between fronts leads to the following type of evolution 
equation for the length of a pulse: 

^ = fci(/.-/i,)-A:4e-^/f + |. (24) 

The coefficients ki and are positive and contain only the real part of the original co- 
efficients. The ffist term is similar to equation ( |16|) where fi^ is the value of the control 
parameter at which a single, isolated front is stationary. Due to the dispersive terms, flc 
includes a correction compared to fic- The second term arises as a result of the overlap of 
the fronts in the convective amplitude. The term involving k^ contains the imaginary parts 
of the coefficients, so that this term represents the interaction due to dispersion. When 
^5 > 0, dispersion provides a repulsive interaction and can lead to the existence of stable 
pulses. When n < fic two pulse solutions exist with the longer one being stable. 

We now consider the combined effects of forcing and dispersion. Of particular interest 
is the question whether periodic forcing can stabilize or destabilize pulses obtained in the 
purely dispersive regime. The task of deriving front equations including both dispersion 
and forcing seems formidable. It is suggested and validated by the numerical simulations 
below that the relevant aspects of both features may be modeled by simply adding the two 
contributions. This leads to an equation of the following form for pulses of length L: 

^ = - + k,v\l - e-t^) - k,v'{l - e~^s^ f - k,e-^'^ + ^. (25) 

If both dispersion and forcing independently result in a stable pulse, their combined 
contribution merely enhances the stability of the pulse. The case of competing contributions 
is more interesting. Figure is a plot of the control parameter /i as a function of steady- 
state pulse lengths L from equation (|25|). The dashed curve is for dispersion in the absence 
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of forcing (z/ = 0), whereas the sohd and dotted curves indicate the addition of forcing. 
For a given value of the control parameter indicated by the thin horizontal line stable 
equilibrium solutions (positive slope) are given by solid circles while open circles indicate 
unstable solutions (negative slope). Plotted are both ) when weak dispersion leads 

to a stable pulse, > 0, and b) when it does not, < 0. If forcing does not stabilize a 
pulse {k2/2k^ > 1), then increasing forcing pushes the stable branch of solutions downward. 
There is a competition between the two interactions and four solution branches may exist 
(Figure [Tl| a inset). Eventually, for sufficiently large forcing z/ the stable pulse disappears 
at infinity, leaving only the unstable pulse solution. Similarly, when dispersion does not 
stabihze the pulse, increasing forcing with k2/2k^ < 1 can eventually lead to a stable pulse 
branch (Figure |ll]b). 

For large L, equation ( p5D is dominated by the dispersive interaction, which decays only 
like l/L, suggesting that the original, dispersive behavior is recovered for large lengths L. If 
k2/2k3 > 1 the competition can lead to two additional pulses, one unstable and one stable, 
for intermediate values of forcing u (cf. inset of Fig. p!T|a). For larger forcing, only one 
unstable pulse and one long, stable pulse remain. Note that for fixed fi and increasing u the 
stable pulse will still disappear at infinity since the branch is being pushed downward. If 
the forcing generates a stable pulse, the 1/L-behavior of the dispersive interaction suggests 



the existence of a long, third, unstable pulse (Fig. |TT]b). Note that according to the analysis 



of the dispersive interaction in |2^, which was confirmed numerically in ||3T[|, the power- 
law behavior continues only up to a maximal length Lmax, beyond which the dispersive 
interaction decays rapidly. The maximal length L^ax decreases with increasing dispersion. 
If forcing dominates the interaction for lengths up to L^ax, then the long-length dispersively 



dominated behavior is no longer expected for either case depicted in Figure |TT|. But it should 
be present for weaker forcing or shorter interaction length. 



Figures |T2| show numerical results that confirm the expected pulse behavior based on 
equation (|25|). Plotted are the numerically obtained pulse lengths as a function of the control 
parameter /i. Stable pulse solutions are indicated by solid symbols and unstable solutions by 
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open symbols. Figures [12|a,b show that increased forcing may destroy a dispersively stable 
pulse. For forcing strength v = 0.2, the forcing is not sufficient to eliminate the stable 
pulse and both solution branches are still present. The inset shows an intermediate value 



of u = 0.4 where four solution branches exist and Figure IT^ shows the pulse lengths for 



stronger forcing z/ = 0.5, where only the unstable branch remains. Figures |T2|c,d depict the 
creation of a stable pulse branch by increasing forcing strength. When u = 0.1 forcing leads 
to a stable branch, but for longer lengths the dispersion dominates and a third, unstable 
pulse exists. However, for u = 0.2 the forcing succeeds in dominating the interaction beyond 
Lmax SO that the branch remains stable. 

VI. CONCLUSIONS 

In this paper we have investigated the effect of external resonant forcing on the interac- 
tion of fronts connecting the stable basic with a stable nonlinear traveling-wave state that 
arises in a subcritical bifurcation. Localized structures were described analytically as bound 
pairs of fronts. The temporal forcing excites the oppositely traveling wave and provides 
an additional mode that is sufficient to localize structures. Since the forcing constitutes 
an externally controlled parameter, pulses of tunable length can be obtained via this new 
localization mechanism. Forcing stabilizes pulses through either a repulsive or an attractive 
interaction between the fronts depending on the pulse length. This is in contrast to other lo- 



calization mechanisms in which stable pulses arise only for a repulsive interaction p7| , |28| , |30 
Multiple pulses with different, but fixed, lengths and holes are also obtained. In addition, 
the combined effect of temporal forcing and dispersion has been investigated. With the 
inclusion of weak dispersion, the interaction between fronts can be described qualitatively 
by a single equation combining the two interaction terms (pSj). It was found that in the 
dispersive regime forcing can lead to the creation of new pulses or the destruction of pulses 
depending on system parameters. The competition between the two interactions determines 
the number and stability of pulses observed. 
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In the regimes investigated here no complex dynamics of the individual pulses have 
been found. It is known, however, that both dispersively stabilized pulses as well as pulses 
stabilized by an advected mode can undergo transitions to chaotic dynamics (e.g. [^9|-|51|). 

Localized traveling waves in the absence of resonant forcing have been observed experi- 
mentally in particular in binary-mixture convection and electroconvection in nematic liquid 
crystals. Depending on parameters, the localization of pulses in binary-mixture convection 
is understood to be due to dispersion and to the coupling to a slowly decaying, advected 
concentration mode. Since the advection depends on the direction of propagation of the 
wave, an interesting question is how the counterpropagating wave excited by a resonant 
forcing will affect the pulses that are stabilized by the concentration mode and how the two 
localization mechanisms interact. The origin of the localization of the worms in electrocon- 
vection is still being investigated. A Ginzburg-Landau model that includes a coupling to 
a weakly damped mode similar to that in the case of binary-mixture convection has been 
proposed By probing the response of the worms to temporal forcing, insight may be 
gained into the relevance of an advected mode in the worms. 

We would like to acknowledge discussions with J. Vifials. This work has been supported 
by the Engineering Research Program of the Office of Basic Engineering Science at the De- 
partment of Energy under grant DE-FG02-92ER14303 and the National Science Foundation 
under grant DMS-9804673. 
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FIGURES 




L 

FIG. 1. Sketch of a pulse traveling to the right. A is the ampUtude of the right-traveling wave. 
Temporal forcing excites the left-traveling wave B, which grows spatially to the left reaching a 
maximum at the trailing front, behind which it decays exponentially. 
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FIG. 2. Sketch of /(A) = A;i(/x - fi^) + k2v'^J^ - k^v'^K^ for fcg > 0. Stable (unstable) fixed 
points are indicated by solid (open) circles. Varying forcing changes the steepness of the parabola 
and the height of the maximum located at A = k2/2k^. Increasing forcing v decreases the stable 
pulse length for (a) /i > /ic and increases the stable pulse length for (b) ^ < fic- 
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FIG. 3. Dependence of pulse length L on the forcing strength v. The symbols indicate the 



numerical results and the curves are obtained from equation (|l9|). The dotted line corresponds 
to = /ic- The curves above (below) this line correspond to values of ^ > /ic < /"c)- The 
dashed curves indicate unstable branches. The group velocity s = 20 and 1.7 in the top and 
bottom figures, respectively, /i = —1.238 (solid squares), —1.2475 (open squares), —1.248075 = 
(circles), —1.2485 (open triangles), and —1.250 (solid triangles). All other parameters are as in 
Figure 13. 
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FIG. 4. Numerically obtained long and short stable pulse solutions, with fi > fic and 
/i < He respectively. The inset zooms in on the amplitude B only. The parameters are 
1/ = 0.2534,5 = 1.7, d = 0.05, c = 2.58, p = 1.0, g = 1.4, r = 4.0, and fi = -1.244 {^l = -1.250) for 
the long (short) pulse. For these parameters /ic = —1.248075. 
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FIG. 5. (a) Pulse velocity as a function of /i for forcing = 0.3, 0.5, and 0.7. The line shown is 
obtained from the analytical result Vpuise = ki{n — ij-c)/"^ + s. (b) Velocity dependence on forcing 
z^. The solid and open symbols indicate the numerically obtained velocities for the trailing and 
leading fronts, respectively. The curves give the analytical results. Squares are for = —1.250 
and circles for fj, = —1.238. The small symbols show the corresponding numerically obtained pulse 
velocities. All other parameters are as in Figure 0. 
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FIG. 6. Amplitude B of pulse solutions when s = 1.25, 5.0, 20.0 indicated by the thick dotted, 
dashed, and solid lines, respectively. The amplitudes A are indicated by thin lines and the leading 
front in A is indistinguishable between the different velocities, /x = —1.2485, v = 0.29 and all other 
parameters as in Figure 0. 
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FIG. 7. (a) Scaled pulse length and (b) relative error of pulse velocity versus forcing u for 
s = 1.25, 5.0, 20.0. fi = —1.2485 and all other parameters as in Figure 0. 



22 



1.2 

1.0 

S 0.8 
ft 

S 0.5 

< 

0.2 



0.0 
35.0 



■ ' 










/ 

A 


B 


w 





45.0 



55.0 



65.0 




10.0 
Hole Length 



20.0 



FIG. 8. (a) Numerical hole solution for c = 1.8, fi = —0.609, and = 0.3. (b) Control parameter 
/i as a function of hole length for various values of u. For c = 1.8, ^3 is in the regime where holes 
exist for both /i > /ic and fj, < fic- Solid (open) symbols refer to stable (unstable) solutions. All 
other parameters are as in Figure . 
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FIG. 10. Two-pulse solution (thick lines) superimposed over a single-pulse solution (thin lines) 
for jj, = —1.238, u = 0.5, and all other parameters as in Figure ^ 
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FIG. 11. Control parameter fi vs. steady-state pulse lengths of equation (p5| ) for dispersion 
only, 1^ = (dashed curve), and added forcing, = 0.5 (solid curve). The solid circles indicate 
stable solutions (positive slope) and the open circles unstable ones (negative slope), (a) For > 
and k2/2k3 > 1 forcing can eventually destroy a stable pulse, d = 0.01 + O.Oi,r = 0.0. The inset 
shows the four branches possible for an intermediate value of forcing ly = 0.425. (b) For k^ < 
and k2/2k3 < 1 increased forcing can lead to stable pulse solutions, d = 0.01 — 0.01i,r = 4.0. All 
other parameters are the same: s = 1.7, c = 2.45 + 0.2i,p = 1.0, g = 1.4. 
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FIG. 12. The control parameter /x vs. the pulse length L. Stable (unstable) solutions indicated 
by solid (open) symbols. Top row (a)-(b) Forcing destroys a dispersively stable pulse {d = 0.01+O.Oz 
and r = 0.0). Forcing strength u = 0.2 and u = 0.4 (inset) in (a), and u = 0.5 in (b). Bottom row 
(c)-(d) Forcing creates a stable branch (d = 0.01 — 0.005z and r = 4.0). Forcing strength u = 0.1 
in (c) u = 0.2 in (d). All other parameters are the same: s = 1.7, c = 2.45 + 0.2i,p = 1.0, g = 1.4. 
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